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ABSTRACT 

We report the first study of Fourier- frequency-dependent coherence and phase/time 
lags at optical wavelengths of cataclysmic variables (MV Lyr and LU Cam) display- 
ing typical flickering variability in white light. Observations were performed on the 
William Herschel Telescope using ULTRACAM. Lightcurves for both systems have 
been obtained with the SDSS filters u' , g' and r' simultaneously with cadences be- 
tween pa 0.5 — 2 seconds, and allow us to probe temporal frequencies between w 10~ 3 
Hz and » 1 Hz. We find high levels of coherence between the u', g' and r' lightcurves 
up to at least « 1CP 2 Hz. Furthermore we detect red/negative lags where the redder 
bands lag the bluer ones at the lowest observed frequencies. For MV Lyr time lags up 
to ss 3 seconds are observed, whilst LU Cam displays larger time lags of ~ 10 seconds. 
Mechanisms which seek to explain red/negative lags observed in X-ray binaries and 
Active Galactic Nuclei involve reflection of photons generated close to the compact 
object onto the surface layers of the accretion disk, where the lag delay is simply the 
light-travel time from the emitting source to the reflecting accretion disk area. Al- 
though this could be a viable explanation for the lags observed in MV Lyr, the lags 
observed in LU Cam are too large to be explained by reflection from the disk and/or 
the donor star. We suggest reprocessing on the thermal timescale of boundary layer 
photons onto the accretion disk as a possible mechanism to explain the lags observed 
in accreting white dwarfs, or reverse (inside-out) shocks within the disk travelling 
through cooler disk regions as they move outwards. 

Key words: accretion, accretion discs - binaries: close - stars: individual: MV Lyr 
and LU Cam - novae, cataclysmic variables - X-rays: binaries 



1 INTRODUCTION 

Cataclysmic variables (CVs) are close interacting binary sys- 
tems where a late-type star transfers material to a white 
dwarf (WD) companion via Roche lobe overflow. With a 
system orbital period ranging from hours to minutes, the 
transferred material from the secondary star forms an ac- 
cretion disk surrounding the WD. As angular momentum is 
transported outwards in the disk, material will approach the 
inner-most regions close to the WD in the absence of strong 
magnetic fields, and eventually accrete onto the compact 
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object. X-ray binaries (XRBs) are also compact interacting 
binaries which are similar to CVs in many ways, but where 
the accreting compact object is either a black hole (BH) or 
a neutron star (NS). Both CVs and XRBs, as well as Active 
Galactic Nuclei (AGN; accreting extragalactic supermassive 
black holes) have been shown to display aperiodic variability 
on a broad range of timescales, with often very strong frac- 
tional rms amplitude. XRBs have shown variability ranging 
from milliseconds to hours, whilst for CVs this ranges from 
seconds to days. This difference can be mainly attributed 
to the fact that the innermost edges of the accretion disks 
in CVs sit at a few thousand gravitational radii (r g ), whilst 
for XRBs it can reach down to just a few r g . The fact that 
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material can get deeper within the gravitational potential of 
XRBs, as compared to CVs, also explains why they are more 
luminous and emit predominantly in X-rays, compared to 
CVs, which emit predominantly at optical/UV wavelengths. 

Aperiodic broad-band variability (also referred to as 
flickering) has been extensively studied in X-rays for XRBs 
ov er several decades in temporal frequency (see for exam- 
plelTerrelj|l972l;lvan der Klisill995l ; iBelloni et alhoool . |2002| ; 
iHoman et al] l200ll ). As CVs emit mostly at optical/UV 
wavelengths, timing studies of these objects had to rely 
on optical observing campaigns from Earth, which are in- 
evitably hindered by large interruptions, poor cadence, and 
in many cases poor signal-to-noise ratios. Furthermore, the 
key timescales to probe in CVs are much longer than in 
XRBs, requiring long, uninterrupted observations. Recently 
CV timing studies have been facilitated thanks to the advent 
of th e Kepler satellite |jenkins et al] l20ld ; ICilliland et al] 
2010), which is able to provide long, uninterrupted and high 
precision lightcurves in the optical light from space. Thanks 
to these capabilities it is now possible to probe over four or- 
ders of magnitude in temporal frequency in CVs. One impor- 
tant discovery in this respect is that the flickering properties 
of CVs in the optical are very similar, at least phenomenolog- 
ically, to those observed in X - rays for XRBs (jScaringi et al] 
l2012bl lal; IWarner et al] 120031 ; iMauche! l2002h . More specifi- 
cally, the discovery of linear rms-flux relations in the CV 
MV Lyr (as well a s in XRBs and AGN lUttlev fc McHardvl 
120011 ; lUttlev et alj|2005h seems to rule out simple additive 
processes as the source of flicker noise (e.g. superposition of 
many independent "shots"), and instead strongly favoured 
multiplicative processes (e.g. mass-transfer variations travel- 
ling from the outer to inner disc for the latter) as the source 
of variability. 

Additional similarities between the broad-band timing 
properties of XRBs and AGN are observed in X-rays when 
analysing simultaneous lightcurves obtained in two energy 
bands. Frequency-dependent phase/time lags are detected 
for both XRBs and AGN over a wide range of frequencies 
dFabian et afll2009l; lNowakll2003l; IVaughan fc Nowaklll997l; 
De Marco et al] I2OI2I ; lUttlev et al] l201ll ; ICassatella et al] 
2012bl a|), where X-ray hard/blue photons lag X-ray soft/red 
ones at low frequencies (by ~ hours for AGN and by w 
seconds for XRBs), and where the opposite is observed at 
higher frequencies (by < 1 hour for AGN and by < 1 second 
for XRBs). Additionally, in some XRBs, the X-ray phase 
lags are observed to change sig n as a function of inten- 
sity,hardness and/or frequencies (|Reig et alj|200cl ). The rea- 
son for these Fourier-dependent phase/time lags is still de- 
bated, but viable scenarios exist to explain the observed phe- 
nomena. For example, the fluctuating acc retion disk model 
l|Lvubarskiilll997l ; lArevalo fc Uttlevll2006l 1 assumes fluctua- 
tions in the mass-transfer rate within the disk to cause the 
observed variability over a wide range of frequencies, with 
long timescale variability being produced further out in the 
accretion disk as compared to the fast timescale variabil- 
ity produced further in. As the long timescale variations 
propagate inwards, they couple to the fast timescale vari- 
ations, which also explains the observed rms-flux relations 
(lUttlev fc McHardvl [200ll ; IScaringi et afll2012bl : iHeil et al] 
120121 ). A consequence of this model is that we should ob- 
serve blue/hard photons lagging redder/softer ones at the 
longest observable frequencies as a consequence of fluctua- 



tions propagating inwards and thus passing through differ- 
ent (and hotter) emitting regions in the disk. On the other 
hand, the soft/negative lag (where soft photons lag hard 
ones) should be observed at the highest frequencies, and are 
explained by reprocessing of X-ray photons produced close 
to the compact object (by the Comptonised component ob- 
served as a power law in hard X-ray spectra) onto the accre- 
tion disk, either as a thermal blackbody (e.g. X-ray heating 
from the disk) or, in the case of AGN, as an additional soft 
photoionised reflection component. In this case the size of 
the negative lags would provide an indication of the size of 
the reprocessing region in the disk: hard/blue photons are 
seen first as they are reprocessed closer in the accretion disk 
than the soft /red ones, and the time lag would simply be 
the light travel time from the central compact object to the 
reprocessing region plus a small reprocessing time. 

A comprehensive analysis of coherence, phase and time 
lags at optical wavelengths has never been performed, and 
because of this no such study has ever been attempted for 
CVfl In this paper we present the first analysis with data 
obtained on the 4.2 meter William Herschel Telescope, using 
ULTRACAM, to study Fourier-dependent coherence, phase 
and time lags over three orders of magnitude in temporal 
frequency for two CVs, LU Cam and MV Lyr. This anal- 
ysis provides useful insight into the broad-band variability 
behaviour of CVs, and reveals further similarities to the be- 
haviour observed in X-rays for XRBs and AGN. 

MV Lyr is classified as being a VY Scl novalike sys- 
tem, spending most of its time in a high state (V~ 12- 
13), but occasionally (every few years) undergoing short- 
duration (weeks t o months) drops in brightness (V«16-18, 
iHoard et alj|2004T ). The reason for these sudden drops in lu- 
minosity is not ccd .. rm Lag lear, but suggestions involving 
star spots on t he secondary inhibitin g mass transfer have 
been proposed (|Livio fc Pringlejll994h . An orbital period of 
3.19 hours has been determined for the system, as well a s 
a low inclinat i on of i « 11° - 13° jSkillman et alj Il995h . 
IScaringi et al] (2012b a) have already studied the single- 
band variability behaviour using Kepler data, and here we 
concentrate on the multi-band behaviour of this system at 
higher frequencies. Much less is known for LU Cam apart 
from its orbi tal period of 3.6 h ours inferred from optical 
spectroscopy l|Sheets et alj|2007l) . We decided to observe LU 
Cam with ULTRACAM to study in more detail the vari- 
ability behaviour after one of us (PJG) noted it strong vari- 
ability on timescal es of weeks in the data of the P alomar 
Transient Factory l|Rau et al]|2009l ; lLaw et al]|2009h . 

In section [2] we will describe our ULTRACAM observa- 
tions and data reduction procedures to obtain the coherence 
and phase/time lags. Section [3] will provide the main results 
and comment on the observed Fourier- dependent features. 
Our discussions and comparison to similar phenomena ob- 
served in XRBs and AGN are presented in section [4] whilst 
our conclusions are drawn in section [5] 



1 Although instruments such as ULTRACAM allow for fast, 
multi-colour and simultaneous photometry, and has been oper- 
ational for several years. 
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2 OBSERVATIONS AND DATA REDUCTION 



2.1 Coherence function 



Both LU Cam and MV Lyr lightcurves presented in this pa - 
per were obtained with ULTRACAM jDhillon et afll2007l l 
at the Cassegrain focus of the 4.2-m William Herschel Tele- 
scope (WHT) on La Palma, Spain. ULTRACAM is a CCD 
camera designed to provide imaging photometry at high 
temporal resolution in three different filters simultaneously. 
The instrument provides a 5 arcminute field on its three 
1024 x 1024 E2V 47-20 CCDs (i.e. 0.3 arcseconds/pixel) on 
the WHT. Incident light is first collimated and then split 
into three different beams using a pair of dichroic beam- 
splitters. For all observations presented here, one beam was 
dedicated to the SDSS u' (3543A) filter, another to the SDSS 
g (4770A) filter and the third to the SDSS r' (623lA) fil- 
ter. Because ULTRACAM employs frame-transfer chips, the 
dead-time between exposures is negligible (« 0.02 seconds). 
The settings used for each observation can be found in Table 
[1] including exposure times and total number of frames. No 
binning was used and the read-out speed was always set to 
slow as it provides the lowest detector readout noise (~ 3 
electrons). 

All data were reduced using the ULTRACAM pipeline 
software. All frames were first debiased and then flat-fielded, 
the latter using the median of twilight sky frames taken 
with the telescope spiralling during twilight before obser- 
vations. For MV Lyr, 2 bright reference stars were used to 
perform differential photometry, whilst 4 bright comparison 
stars were used for LU Cam. Furthermore, LU Cam was ob- 
served in the u'-band in co-add mode. For this mode, one 
tt'-band exposure is performed for every two g' and r' expo- 
sures to increase signal-to-noise, but nevertheless all three 
exposures remain synced. To place all filters on the same 
time stamp, we averaged pairs of the g' and r' exposures. 
Portions of the extracted lightcurves for MV Lyr and LU 
Cam in the three filters are shown in Fig[T] 

For each observing run (two for MV Lyr and one for 
LU Cam) we obtained time averaged power spectral den- 
sities (PSDs) by splitting each observing run (each lasting 
~ 6 hours) into 6 non-overlapping segments. For each seg- 
ment (and for each band) we computed the Fast-Fourier 
transform (FFT) a nd ap plied the rms normalisation of 
iBelloni fc HasingeJ (| 19901 ) so that the square root of the 
integrated PSD power over a specific frequency range yields 
the root-mean-square (rms) variability. Fig. [2] shows the 
computed PSDs (in power x frequency and with Poisson 
noise subtracted) for all objects and in all bands. All are 
roughly described with a broken powerlaw up to ~ 10 _1 
Hz, after which Poisson noise dominates the PSDs. The PSD 
br eak at sa 10~ 3 Hz for M V Lyr is consistent with the results 
of IScaringi et~afl |2012bl lah, and discussed therein. 

Additionally, we carried out analysis of cross-spectral 
Fourier statistics. These are the coherence, phase and 
time lags, which have been applied to X-ray timing data 
of XRBs and AGN in the past, and are described in 
IVaughan fe Nowakl (|l997h and lNowakl |2003l ). Because this 
is the first time such methods are applied to optical data, 
we here review all the steps involved again, with particular 
emphasis on optical observations of CVs. 



The coherence function, j 2 , is a Fourier- frequency- 
dependent measure of the degree of linear correlation be- 
tween two simultaneously observed lightcurves in two differ- 
ent filters (or energy bands as used in X-rays). Specifically, 
it provides the fraction of the mean-squared variability at 
a specific frequency of one lightcurve which can be directly 
attributed to the other by a linear transformation (e.g. a 
simple time shift and flux re-scaling). As a simple intuitive 
example, two simultaneously observed lighturves in two fil- 
ters which display a constant time shift relating one to the 
other will display perfect coherence as well as a constant 
time lag, at all frequencies. 

Let us consider two simultaneously observed lightcurves 
in two different filter s r(t) and u(t) . We here use a similar no- 
tation to that used in iNowakl (|2003h . so that r(t) is the longer 
wavelength band and u(t) the short wavelength band (s(t) 
and h(t) respectively in lNowak1l2003h . To compute the coher- 
ence (and indeed any other higher-order Fourier statistic) we 
require an ensemble of independent measurements for each 
band in order to reduce statistical noise. Thus both r(t) and 
u(t) are split into M independent segments, where M — 6 
in our case for both LU Cam and MV Lyr, so that each seg- 
ment is about 1 hour long. For each lightcurve segment, 
ri=i->-M(t) and Ui=i^-M (t), we compute the Fast Fourier 
transforms (FFT) X r< i(f) and X Ui i(f) respectively, together 
with the corresponding power spectra P r ,i(f) = \X ri i\ 2 and 
Pu,i(f) = \X U)i \ 2 and cross-spectra C»(/) = X*^(f)X u ,i(f) 
(where * denotes the complex-conjugate). 

For each FFT the Poisson (white) noise levels can be 



calculated as |iV r , 



and \N U 



for both 



bands (using the lBelloni fc Hasingerl |l99(j) normalisation). 
Because in our case the errors obtained from the differen- 
tial photometry will not be Poissonian (and could possibly 
be correlated between the different filter lightcurves due to 
systematic seeing effects in the comparison stars) we fit each 
independent PSD with a powerlaw plus a constant to esti- 



mate the mean white noise levels \N r 



and \Nu,i 



The 



intrinsic power for each segment in the r lightcurve is then 
defined as |-Ri| 2 = P r ,i — \N ri i\ 2 , and similarly for the u 
lightcurve as \U l \ 2 = P Ull - \N Uri \ 2 . 

From this, we first define the raw coherence (i.e. the co- 
herence without taking into account white noise) at a spe- 
cific frequency as 



2 

'Yraw 



(/) 



l«C(/)»| 



l«C(/))>| 



«|X r (/)P»«|X u (/)p» «P,.(/)»«P u (/)» 

(1) 

where the double angled brackets denote averages over ad- 
jacent frequency bins considered (n/ s in total) as well as av- 
erages between the M independent segments. This notation 
will be used throughout the paper. The frequency binning 
in our case is logarithmic, with the constrain that each bin 
requires at least 5 datapoints, thus n/ s ^ 5 for all bins. The 
associated error is then defined as 



^^iraw (y ) 



2 l laJJW ~ 1ra W (f)} 
\lraw(f)\ 



(2) 



where m = M x rif s , the number of measurements used in 
a specific frequency bin. 

The intrinsic coherence, ■y 2 , can be estimated by cor- 
recting each term in Eq. [1] for counting noise. The term 
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Figure 1. ULTRACAM lightcurves obtained in the u' (bottom- blue), g' (middle-green) and r' (top-red) filters for the CVs MV Lyr 
(top two panels) and LU Cam (bottom panel). The lightcurves have all been normalised by their respective mean count rate. The r' and 
u' lightcurves have also been shifted by 0.2 and —0.2 respectively for clarity. Both systems clearly display variability on a wide range of 
frequencies. Colour figure is available online. 



Object Run Dates Exposure/frame (s) Frames 

MV Lyr I 31 August 2012 0.844 27735 

MV Lyr 2 1 September 2012 0.526 41670 

LU Cam 3 10 January 2012 2.276 18502 

Table 1. ULTRACAM journal of observations during the observing runs. 



|((C(/)})| 2 contains a positive bias (offset) caused by Pois- 
son noise, which can be estimated as 

n 2 = ({(\R\ 2 )){(\N u \ 2 )) + {(\U\ 2 ))({\N r \ 2 )) (3) 

+<(W 2 »«|AT u | 2 »)/m. 
With this, we can define the intrinsic coherence by 

|((C(/)))| 2 -n 2 



YU) = 



{{\RU)\ 2 )) {{\u{f )?))'■ 



(4) 



where the denominator now includes only intrinsic power. 
The associated error is then defined as 



2n m 



<57 2 (/)= 7 2 (/) 



, «l^-| 4 )) , ((\N U \ 4 )) m5^ av 



(K{C)) |2_ n2) 2 {{lm) {{m 



2.2 Fourier phase/time lags 



m 

1/2 



The Fourier phase/time lags are also computed between 
two simultaneous lightcurves obtained in different wave- 
length ranges, and like the coh erence function, are re- 
lated to the cro s s-corr e lation dBendat fc Piersoll [2010; 
IVaughan fc Nowakl 1 19971 ; iNowakl boost ). Similarly to the 
PSDs in Fig[2] and the coherence, we split the lightcurves 
into M — 6 independent segments. The Fourier phase lag is 
defined as the phase of the average cross power spectrum, 
defined as 



*(/) = arg[((C(f)))}, 



with associated errors defined as 

SHI) -- ' 1 



-y 2 

I raw 



2^ 2 

^ ira'u 



(7) 



The Fourier time lag is constructed by simply dividing 
the phase lag by 271"/, so that r(f) — Similarly, the 

error on the time lag is defined as Sr(f) = -f^p- 



3 RESULTS 

(5) Fig.[3]shows the coherence functions, phase and time lags for 
the three observing runs and all three colour combinations 
v! — r', it' — g' and g' — r' . Both LU Cam and MV Lyr show 
high levels of coherence at the lowest observed frequencies in 
all three colour combinations. We note that above w 10 -2 Hz 
the noise levels start to affect our observations, resulting in 
larger error bars in this frequency range. The coherence is a 
fourth-order statistic and is very sensitive to uncertainties i n 
the noise subtraction (jBendat fc Piersolll2010l ; lNo wak 2003). 

The noise levels we infer from the PSD fits may also 
be slightly over/underestimated for our observations as the 
errors are not strictly Poissonian. Specifically, we expect er- 
rors in adjacent filter lightcurves to be correlated to some 
degree due systematic effects in the comparison stars. For 
example, the apparent rise above 1 in the intrinsic coher- 
ence plots of MV Lyr in the r' — g' combination are due 
to an overestimation of the white noise level. Because of 
(6) these uncertainties and possible correlations within errors 
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Figure 2. Rms normalised PSDs (power X frequency with Poisson noise subtracted) for LU Cam (left) and MV Lyr (right, 2 different 
nights) in the u 1 (blue squares), g' (green circles) and r' (red diamonds) filters. Each PSD has been obtained by averaging 6 independent 
PSDs, each obtained from pal hour timeseries. Both obje cts show a flattening o f the variability power below ss 10 — 3 Hz. In the case of 
MV Lyr, this is consistent with the Kepler observations of Scaringi et al. ( 20123H) ■ Also, both objects display higher levels of variability 
with increasing photon energy. Colour figure available online. 



of the different colour combination, the apparently signifi- 
cant features appearing in multiple band selections at high 
frequencies must be treated with caution. Nevertheless we 
are confident our coherence and phase/time lag estimates 
are reliable for frequencies below « 10 -2 Hz where uncer- 
tainties in the noise levels are negligible compared to the 
variability power, as also illustrated by the raw coherence 
plots. 

Both LU Cam and MV Lyr also show similar features 
in the phase/time lags below ~ 10~ 2 Hz. Specifically, we 
detect significant red/negative lags for both systems, where 
low energy photons are seen arriving after the high energy 
ones. In the case of LU Cam negative lags are observed with 
a «10 second delay below 10 -3 Hz, whilst for MV Lyr this 
is much smaller at ~3 seconds. Furthermore in both systems 
the negative lags are most pronounced in the r' — v! com- 
bination, which might be expected as this is the filter com- 
bination with the largest wavelength difference. At higher 
frequencies, the time lags in MV Lyr seem to asymptote to 
0, whilst for LU Cam they seem to swap sign reaching a 
few seconds at the highest observable frequencies. However, 
because the Poisson noise level in this frequency range af- 
fects our observations as described above, the high frequency 
results must be treated with caution. 



4 DISCUSSION 

As this is the first study of this kind in CVs, it is non- 
trivial to predict the Fourier dependent behaviour we have 
observed. We will thus address our results in the context 
of different timescales which could potentially explain the 
observed lags in MV Lyr and LU Cam. 



4.1 Viscous timescale 

If we take the fluctuating accretion disk model, which is fre- 
quently invoked to explain the Fourier time l ags observed 
in the X-ray domain for XRBs and AGN (|Uttlev et al.l 



201ll; ICassatella et al l | 2012bl lal; lUttlev fc McHardvl l200ll ; 
Arevalo k, Uttlevl 12006 ). we would expect to observe 
blue/positive lags in the optical for CVs. This would be 
a consequence of fluctuations propagating inwards through 
the disk on the viscous timescale, and emitting redder pho- 
tons before the blue ones as a consequence of accretion rate 
fluctuations passing through lower and then higher temper- 
ature emitting regions as they move inwards through the 
disk. What we observe in MV Lyr and LU Cam is the oppo- 
site: variations in blue photons are observed preceding the 
variations in the red photons. Because of this we must rule 
out this timescale as a possible explanation for the observed 
lags. 



4.2 Light-travel timescale 

Red/negative lags somewhat similar to those observed 
here have also been observed in X-rays for XRBs and 

AGN (referred to as soft /negative lags). For example, 

|Pe Marco et all l|2012l ) have shown how soft/negative lags 
are observed in X-rays for a sample of AGN, up to hundreds 
of seconds. The main interpretation for these lags is that 
they represent the light travel time to the disc from the vari- 
able continuum source. The continuum source illuminates 
the disc, which in the case of AGN leads to a soft pho- 
toionised reflection spectrum, while in XRBs the absorbed 
photons are also reprocessed and re-emitted as (soft X-ray) 
thermal radiation at the local disc blackbody temperature. 
Thus, in this model, the delay is simply the light- crossing 
time (plus a small reprocessing time, assuming also that the 
time-scale for thermal reprocessing is small, see Section r4.3p . 
which for XRBs will be milliseconds, and for AGN will be on 
the order of tens to hundreds of seconds. This scenario could 
potentially explain the ~ 1 — 3 second lags observed in MV 
Lyr, placing the reflecting/reprocessing region between 0.4- 
0.8J?q (if we take the sy stem parameters from lLinnell et all 
120051 ; iHoard et al.ll2004l ). This range is both consistent with 
the outer-edges of the disk and/or the secondary star. 

Although the lag values are consistent with the expected 
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Figure 3. Coherence and phase/time lags observed with ULTRACAM for LU Cam (left) and MV Lyr (right, 2 different nights) in the 
three combinations r' — g' (red diamonds), g' — u' (blue circles) and r' — u' (black squares). All have been obtained from averaging 6 
independent lightcurve segments. Negative time lags indicate that the emission in the redder bands lags the emission in the bluer bands. 
Colour figure available online. 



binary separation in MV Lyr, this model would require the 
r-band emission to be dominated by the reprocessed light, 
otherwise the lags would be diluted by the variable central 
continuum emission. Furthermore, this interpretation seems 
to be unlikely for disc reprocessing at these large radii, since 
the disk is unli kely to reprocess much of the central contin- 
uum emission (|van Paradiis fc Verbundll984l ; Ide Jong et al.l 
1 19961 ). Additionally, the much larger lags of ~ 10 seconds 
observed in LU Cam rule out the possibility of the lag being 
produced by light-travel times to the outer disc and/or sec- 
ondary star. LU Cam does however leave the possibility that 
the reflecting region would sit at large radii, at ~ 6Rq, out- 
side of the binary orbit, maybe in the form of circumbinary 
disk/torus, although we would have to be viewing the system 
from a very specific angle. If this were the case, we would 
also expect the large scale varying emission to be weak com- 
pared to the direct continuum emission from close to the 
compact object. This would imply that the observed lags 
from the reprocessing region at large distances would be di- 
luted from the direct continuum emission, meaning that the 
intrinsic lags should be even larger than the ones observed 
in LU Cam. We note however that for the possibilities pre- 
sented here further spectral analysis would be required to 
investigate these possibilities. 



4.3 Thermal timescale 

Another, possibly more plausible, scenario could be that the 
accretion disk is reprocessing energy originating from close 
to the compact object, for example in the boundary layer, on 
the local thermal timesc ale. In MV Lyr, the b oundary layer 
can reach « 100, 000 K jGodon fc Sionll201ll ), allowing the 
colder disk (of w 10, 000 K in the outer edges) to efficiently 
absorb such photons (as opposed to reflect them through 
Compton scattering as is the case for XRBs). In this scenario 
one can imagine variable radiation influencing the accretion 
disk, where the observed red/negative lags are due to the 
light-crossing time from the source to the disk, but also due 
to the thermal reprocessing time of photons within the disk 
before they are re-emitted and observed. In this case the 
relevant timescale to consider is the thermal timescale 



where SI is the Keplerian frequency at a specific disk radius, 
and a the disk viscosity paramet er, also at a specific disk 
radius (|Shakura fc Sunvaevlll973l ). With an accretion disk 
reaching the WD surface of « O.Oli?©, the corresponding 
thermal timescale at the inner-most edges of the disk is on 
the order of tens of seconds, depending on a. Our LU Cam 
result could potentially be explained through this process if 
a ~ 0.7, which is considered to be high for accretion disks in 
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CV systems. If on the other hand the reprocessing of photons 
were to occur at larger radii (say close to the outer edges 
of the disk at m 0.5R@), then the thermal timescale would 
be much longer. In both scenarios however we would only 
expect the surface layer of the disk to reprocess photons, 
greatly reducing the timescale inferred from Eq. [5] and pos- 
sibly explaining the observed lags in MV Lyr and LU Cam. 
We note that the above scenario is a very general example, 
and we can imagine magnetic fields (instead of photon radi- 
ation) affecting the disk, which reprocesses the energy, and 
emitting photons after some delay, but again comparable to 
the thermal timescale. The larger red/negative lags observed 
in AGN are also explained by a similar process, where the 
inner-disk edge lies at a few AU. The main difference is that 
for AGN (and XRBs) the X-ray photons originating close 
to the compact object are re-emitted nearly instantaneously 
as they photoionise the surface layers of the disk. In CVs 
on the other hand the optical/UV emission close to the WD 
cannot photoionise the disk, but instead heats it up, and 
photons from the disk are then re-emitted on the thermal 
timescale. 

One last possible explanation for the observed lags in 
LU Cam and MV Lyr is th at of reverse (inside-o ut) shocks 
within the accretion disk (|Krauland et al.l l2013h . possibly 
originating from the WD boundary layer accretion disk. 
These waves would then transport energy outwards in the 
disk. In this scenario the lightcurve variations in the u band 
will be observed before those in the r band as a consequence 
of the shock(s) moving outwards in the disk and passing 
through the hotter inner-edges and later throu gh the cooler 
outer -edges. The thermal timescale of the disk ()Frank et al.l 
120021 ) would also be the relevant timescale for the propa- 
gation of the Shockwave, potentially in line with our obser- 
vations. This interpretation, as well as all others presented 
here, require thorough modelling in order to explain the lag 
phenomena observed in MV Lyr and LU Cam. Such mod- 
elling which however beyond the scope of this paper. 



5 CONCLUSION 

We have presented the first analysis of coherence, phase and 
time lags for CVs at optical wavelengths using ULTRACAM, 
mounted on the William Herschel Telescope. Our analysis is 
based on two objects, namely MV Lyr and LU Cam, which 
both show significant red/negative lags at frequencies be- 
low « 10 -2 , where blue photons are observed before the red 
ones. For MV Lyr this lag is observed to a maximum of « 3 
seconds, whilst for LU Cam the lag is much larger at w 10 
seconds. Furthermore we established that both objects show 
high levels of coherence over a wide range of frequencies up 
to at least « 10~ 2 Hz. Mechanisms to explain the observed 
phenomena have been proposed, involving reprocessing of 
boundary layer photons (or other sources of energy close to 
the compact object) from the inner-edge of the accretion 
disk on the local thermal timescale, or reverse (inside-out) 
shocks travelling outwards in the disk. A thorough study 
of the implications of the results presented here, together 
with detailed modelling of the accretion disk reprocessing is 
beyond the scope of this Paper. We note however that any 
models trying to reproduce the observed Fourier-dependent 
behaviour in XRBs and AGN will, in future, need to ad- 



dress the features presented here as well, in order to provide 
a self-consistent model explaining the variability behaviour 
through all wavelengths (including optical) and throughout 
the different compact accretors (including white dwarfs). 
Here we have only presented results on two CVs, and noted 
that both display significant lags at optical wavelengths. We 
might also expect that most flickering CVs will display simi- 
lar lags, which will provide further grounds on which to test 
models seeking to explain the variability properties observed 
throughout the different types of compact accreting objects. 



ACKNOWLEDGEMENTS 

This research has made use of NASA's Astrophysics Data 
System Bibliographic Services. S.S. acknowledges funding 
from the FWO Pegasus Marie Curie Fellowship program, as 
well as funding from NWO project 600.065. 140.08N306 to 
P.J.G. S.S. also wishes to thank Pablo Cassatella for useful 
and insightful discussions. E.K. acknowledges funding form 
the NWO Vidi grant 639.042.218. 

REFERENCES 

Arevalo P., Uttley P., 2006, MNRAS, 367, 801 

Belloni T., Hasinger G., 1990, A&A, 230, 103 

Belloni T., Klein-Wort M., Mendez M., van der Klis M., 

van Paradijs J., 2000, A&A, 355, 271 
Belloni T., Psaltis D., van der Klis M., 2002, ApJ, 572, 392 
Bendat J., Piersol A., 2010, Random Data: Analysis and 

Measurement Procedures, Wiley Series in Probability and 

Statistics. Wiley 
Cassatella P., Uttley P., Maccarone T. J., 2012a, ArXiv 

e-prints 

Cassatella P., Uttley P., Wilms J., Poutanen J., 2012b, MN- 
RAS, 422, 2407 

de Jong J. A., van Paradijs J., Augusteijn T., 1996, A&A, 
314, 484 

De Marco B., Ponti G., Cappi M., Dadina M., Uttley P., 
Cackett E. M., Fabian A. C, Miniutti G., 2012, ArXiv 
e-prints 

Dhillon V. S., Marsh T. R., Stevenson M. J., et al., 2007, 

MNRAS, 378, 825 
Fabian A. C, Zoghbi A., Ross R. R., Uttley P., Gallo L. C, 

Brandt W. N., Blustin A. J., Boiler T., Caballero-Garcia 

M. D., Larsson J., Miller J. M., Miniutti G., Ponti G., 

Reis R. C, Reynolds C. S., Tanaka Y., Young A. J., 2009, 

Nature, 459, 540 
Frank J., King A., Raine D. J., 2002, Accretion Power in 

Astrophysics: Third Edition. Cambridge Press 
Gilliland R. L., Jenkins J. M., Borucki W. J., et al., 2010, 

ApJL, 713, L160 
Godon P., Sion E. M., 2011, PASP, 123, 903 
Heil L. M., Vaughan S., Uttley P., 2012, MNRAS, 422, 2620 
Hoard D. W., Linnell A. P., Szkody P., et al., 2004, ApJ, 

604, 346 

Homan J., Wijnands R., van der Klis M., Belloni T., van 
Paradijs J., Klein-Walt M., Fender R., Mendez M., 2001, 
ApJS, 132, 377 

Jenkins J. M., Caldwell D. A., Chandrasekaran H., et al., 
2010, ApJL, 713, L87 



8 S. Scaringi et al. 



Krauland C. M., Drake R. P., Kuranz C. C, Loupias B., 
Plewa T., Huntington C. M., et al., 2013, ApJL, 762, L2 

Law N. M., Kulkarni S. R., Dekany R. G., Ofek E. O., 
Quimby R. M., Nugent P. E., et al., 2009, PASP, 121, 
1395 

Linnell A. P., Szkody P., Gansicke B., et al., 2005, ApJ, 
624, 923 

Livio M., Pringle J. E., 1994, ApJ, 427, 956 
Lyubarskii Y. E., 1997, MNRAS, 292, 679 
Mauche C. W., 2002, ApJ, 580, 423 

Nowak M. A., 2003, in New Views on Microquasars, 

P. Durouchoux, Y. Fuchs, & J. Rodriguez, ed., p. 11 
Rau A., Kulkarni S. R., Law N. M., Bloom J. S., Ciardi D., 

Djorgovski G. S., et al., 2009, PASP, 121, 1334 
Reig P., Belloni T., van der Klis M., Mendez M., Kylafis 

N. D., Ford E. C, 2000, ApJ, 541, 883 
Scaringi S., Kording E., Uttley P., Groot P. J., Knigge C, 

Still M., et al., 2012a, MNRAS, 427, 3396 
Scaringi S., Kording E., Uttley P., Knigge C, Groot P. J., 

Still M., 2012b, MNRAS, 421, 2854 
Shakura N. I., Sunyaev R. A., 1973, A&A, 24, 337 
Sheets H. A., Thorstensen J. R., Peters C. J., Kapusta 

A. B., Taylor C. J., 2007, PASP, 119, 494 
Skillman D. R., Patterson J., Thorstensen J. R., 1995, 

PASP, 107, 545 
Terrell Jr. N. J., 1972, ApJL, 174, L35 
Uttley P., McHardy I. M., 2001, MNRAS, 323, L26 
Uttley P., McHardy I. M., Vaughan S., 2005, MNRAS, 359, 

345 

Uttley P., Wilkinson T., Cassatella P., Wilms J., 
Pottschmidt K., Hanke M., Bock M., 2011, MNRAS, 414, 
L60 

van der Klis M., 1995, X-ray Binaries, 252 

van Paradijs J., Verbunt F., 1984, in American Institute of 
Physics Conference Series, Vol. 115, American Institute of 
Physics Conference Series, Woosley S. E., ed., pp. 49-62 

Vaughan B. A., Nowak M. A., 1997, ApJL, 474, L43 

Warner B., Woudt P. A., Pretorius M. L., 2003, MNRAS, 
344, 1193 



